UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADP012036 

TITLE:  Approximated  Planes  in  Parallel  Coordinates 
DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 


This  paper  is  part  of  the  following  report: 

TITLE:  International  Conference  on  Curves  and  Surfaces  [4th],  Saint-Malo, 
France,  1-7  July  1999.  Proceedings,  Volume  1.  Curve  and  Surface  Design 

To  order  the  complete  compilation  report,  use:  ADA399461 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADPO 1 20 1 0 thru  ADPO 1 2054 


UNCLASSIFIED 


Approximated  Planes  in  Parallel  Coordinates 


Tanya  Matskewich,  Alfred  Inselberg  and  Michel  Bercovier 


Abstract.  For  the  visualization  of  multivariate  problems,  a multidimen- 
sional system  of  Parallel  coordinates  is  used  which  provides  a one-to-one 
mapping  between  subsets  of  N-space  and  subsets  of  2-space.  A rigorous 
methodology  for  doing  and  seeing  N-dimensional  geometry  emerges  as  well 
as  several  applications.  Here  an  application  to  Error  Tolerancing  involv- 
ing the  visualization  and  characterization  of  “approximate  coplanarity”  is 
presented.  The  exact  description  of  the  neighborhood  of  an  TV-dimensional 
hyperplane  in  a parallel  coordinate  system  is  given. 


§1.  Introduction 

The  parallel  coordinate  system  serves  as  a tool  for  visualization  of  multi- 
dimensional objects  and  multivariate  relations.  It  was  shown  ([1,2])  that  this 
representation  gives  a simple  and  constructive  geometrical  description  for  sub- 
sets of  points  which  are  strictly  coplanar  (i.e.  belong  to  a common  p-flat  in 
TV-dimensional  case).  This  allows  the  visualization  of  coplanar  points  and 
the  existence  of  linear  dependencies  between  variables.  It  leads  to  numerous 
applications  in  different  fields,  and  also  practical  applications  involving  finite 
error  tolerancing.  Here  exact  descriptions  of  approximated  hyperplanes  in 
the  parallel  coordinate  system  are  given,  providing  a methodology  for  their 
visualization. 

In  the  next  section  a brief  review  of  previous  results  is  given.  It  is  followed 
in  Section  [3]  by  an  exact  mathematical  formulation  of  the  problem  in  the 
general  case.  Sections  [4]  and  [5]  contain  some  auxiliary  lemmas  which  make 
the  main  result  more  intuitive,  as  well  as  the  main  result  itself. 

Although  we  have  a complete  and  precise  proof  of  the  main  result  in  the 
general  case  of  “approximated”  p- flats  in  TV-dimensional  space,  lack  of  space 
prevents  us  from  presenting  it  here. 
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Fig.  1.  Three  collinear  points  in  2D  (left)  and  3D  (right). 


§2.  Representation  of  Affine  Subspaces  in  Parallel  Coordinates 

The  parallel  coordinate  system  is  constructed  in  the  following  way:  in  the 
Euclidean  plane  R2  (or  more  precisely  in  the  2D  projective  plane  P2  ) with 
xj/-Cartesian  coordinates,  N copies  of  the  axis  y labeled  X\, ...  ,X n are  placed 
equidistant  (usually  the  distance  between  adjacent  axes  is  taken  as  1)  and 
perpendicular  to  the  x-axis.  They  are  the  axis  of  the  Parallel  Coordinate 
system. 

A point  with  Euclidean  coordinates  (pi, . . ■ ,pn)  is  represented  by  a polyg- 
onal line  with  N vertices  fj  = (i  — l,p,),  one  on  each  axis.  In  this  way,  a 1 1 

correspondence  between  points  in  R.,v  and  planar  polygonal  lines  with  vertices 
on  the  parallel  axes  is  established. 

In  2D,  a point  is  represented  by  a line  (usually  just  the  segment  between 
the  axes  is  shown).  It  can  be  easily  proved  that  the  lines  representing  points 
of  a line  C\X\  + C2X2  = co  (for  Ci  + c 2 ^ 0)  intersect  at  the  point  7fi2  with 
xy-coordinates  ( cf*C2 , ~Clc+C2 ) 1 or  more  generally  at  the  point 


/ d\C\  + d2C2  Cq  \ 

\ ci  + c2  ’ ci  +c2)  ’ 


(1) 


where  d\  and  d.2  are  distances  between  axis  y and  A'i,  X2,  respectively.  (Lines 
with  slope  1 are  mapped  onto  the  ideal  points  of  projective  plane,  but  in  what 
follows  we  will  not  consider  any  “degenerate”  cases).  Hence  a fundamental 
point  <->  line  duality  is  induced  (see  Fig.  1). 

In  3 D,  a line  can  be  fully  described  by  any  pair  from  its  three  projec- 
tions on  coordinate  planes.  Each  such  projection  is  a line  in  2D-space  of  the 
corresponding  coordinates,  and  so  can  be  represented  in  parallel  coordinates 
exactly  as  was  described  above.  Hence,  if  4 — Cq‘j^  - projection 

of  the  line  on  Xi Xj  Euclidean  plane  (i,j  € {1,  2,  3}),  then  it  is  represented  in 
the  parallel  coordinate  system  by  the  point  tr ij  whose  coordinates  may  be  com- 
puted from  equation  (1),  and  can  be  found  geometrically  as  the  intersection 
of  corresponding  lines.  The  three  points  7fi2,  tr23  and  7fi3  are  always  collinear 
as  a consequence  of  Desargue’s  Theorem,  and  any  two  of  them  represent  the 
line  in  parallel  coordinates.  We  denote  by  L the  line  on  which  the  three  points 
lie  (see  Fig.  1). 


Approximated  Planes  in  Parallel  Coordinates. 


259 


Fig.  2.  Randomly  sampled  set  of  coplanar  points  in  3D  and  intersection  of  lines 
L showing  coplanarity. 

Let  us  now  consider  a plane  C1X1  + c2x2  + C3X3  = cq.  Every  pair  of  points 
belonging  to  this  plane  define  a line  L in  parallel  coordinates  (which  can  be 
constructed  geometrically  from  the  representation  of  points  themselves).  It 
can  be  shown  by  direct  computation  that  all  such  lines  intersect  at  the  common 
point  with  coordinates 


7T123  = 


CO 


ci  + 2c2 


ci  + c2  + c3  Cl  + c2  + c3 


(2) 


This  condition  can  be  used  to  characterize  coplanarity.  Note,  that  one  point 
is  not  sufficient  in  order  to  specify  the  plane.  The  solution  is  to  introduce 
an  additional  axis  X[  placed  after  X3  and  at  a unit  distance  from  it,  and 
consider  the  representation  of  points  also  in  axes  (X2,X3,X[).  This  leads  to 

the  additional  point  7f23i'  = ^+g+^)  (see  FiS-  2)- 

This  generalizes  nicely  to  the  iV-dimensional  case,  and  it  can  be  shown 
that  a representation  of  a hyperplane  in  parallel  coordinates  also  can  be  re- 
cursively constructed  by  a simple  geometric  procedure,  using  affine  subspaces 
of  lower  dimensions.  A hyperplane  is  represented  by  N — 1 indexed  points; 
the  “first”  one  has  coordinates 


7r12...iV  = 


c2  + 2c3  + . . . + (N  - l)cjy 
ci  + c2  + . . . 4-  c/v 


co 


ci  + c2  + . . . + c/v 


(3) 


and  the  others  have  very  similar  formulas. 

A p-flat  in  TV-dimensional  case  can  be  described  by  N — p linearly  inde- 
pendent equations,  where  each  of  them  has  the  form  1 cikxit  = co,  and 
so  corresponds  to  a hyperplane  in  axes  (A), , Xi2, . . . , Xip+1).  It  follows  that 
a p-flat  is  represented  by  p(N  — p)  indexed  points.  The  ensuing  discussion  is 
restricted  to  “approximated”  hyperplanes.  It  is  easy  to  show  that  the  general 
case  of  an  “approximate”  p-flat  in  TV-dimensional  space  can  be  reduced  to  the 
study  of  some  “approximate”  hyperplane. 
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Fig.  3.  Sampled  proximate  line  ax  4-  by  = 1,  a £ [a  , a+],  b € [6  , 6+]  in 
orthogonal  and  parallel  coordinates. 


§3.  Exact  Problem  Definition 

We  will  use  the  following  definition  of  “similarity” : 

Definition  1.  Proximate  flats  are  defined  here  as  flats  with  proximate  equa- 
tions. An  approximate  hyperplane  is  defined  as  a set  of  hyperplanes  given  by 
equations 

C\X\  + C2X2  + • • • + CtfXN  = Co,  (4) 

where  every  coefficient  can  vary:  c*  € [c~ , cf] , i = 0, . . . , N. 

Such  a slab  of  hyperplanes  is  extremely  difficult  to  visualize  in  orthogonal 
coordinates,  even  for  3 D.  Another  problem  is  that  even  in  the  2D  case,  line 
neighborhoods  are  unbounded  in  orthogonal  coordinates,  so  neighborhoods  for 
different  lines  always  overlap.  Fig.  3 shows  that  samples  of  proximate  lines  (in 
2D)  form  a cloud  in  the  form  of  a very  simple  and  nice  convex  quadrilateral. 

In  the  TV-dimensional  case,  a hyperplane  is  described  by  N — 1 points 
which  means  that  in  the  approximated  case,  we  will  get  N — 1 “clouds”  of 
indexed  points  in  parallel  coordinates.  In  order  to  make  things  simpler,  we 
will  use  the  following 

Assumption  2.  Free  coefficients  of  equations  are  not  allowed  to  vary,  and 
are  supposed  to  be  identically  equal  to  1. 

(We  have  a complete  analogue  of  the  main  result  for  the  case  when  this  as- 
sumption is  not  applied.) 

Lemma  3.  It  is  sufficient  to  study  only  the  range  of  the  first  indexed  point 
of  a hyperplane.  That  is,  the  following  mathematical  problem  should  be  con- 
sidered: find  the  range  of  the  function  f : IRV  — > ]R2  such  that 


/ (ci , • • • ,Cjv) 


EjLiO'-ite  1 \ 

yN  ’ yN  J 

A—sj  — 1 Z^j-1  CJ  / 


(5) 


when  cj  € [cj  ,c+]  (cj  < cf,  j = 1, . . . , N). 
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Fig.  4.  Possible  place  for  domain  f l depending  on  the  signs  of  CoeffMin  k and 
CoeffMaxf.. 

The  notations  used  are: 

• a:(ci, . . . , cn)  and  y(c\, ...  ,cn)  - the  first  and  the  second  coordinates  of 
/jv(c i, . . . , cjv)  in  1R2  respectively; 

• B = [c^- , cf  ] x . . . x [c“,c^]  - box  in  the  space  of  coefficients; 

• fl  = f(B)  - image  of  B in  the  parallel  coordinate  system. 

We  now  show  why  Lemma  3 holds  in  the  3 D case.  Here  studying  the  range 
of  ^123  (equation  (2))  is  sufficient  because  ir23i<  can  be  rewritten  in  the  form 

= (1,0)  + ( -cl+2c'  , — 1— ■ ) 

\Ci  + c2  + C3  Cl  + c2  + c3  / 

This  implies  that  the  equation  of  if23i<  can  be  obtained  from  the  one  of  7fi23 
by  shift  and  cyclical  change  of  parameters  ci  -+  c2,  c2  — ♦ C3,  C3  — > ci. 

For  the  general  case,  this  “reduction”  lemma  can  be  proved  using  simple 
combinatorics,  and  it  can  be  shown  that  in  order  to  describe  the  range  of  one 
of  the  indexed  points  corresponding  to  a p-flat  in  iV-dimensional  space,  it  is 
always  sufficient  to  consider  /2jy  with  some  coefficients  identically  equal  to 
zero. 


§4.  Some  Notes  on  the  Domain  fl 

To  get  an  intuitive  feel  about  the  structure  of  the  domain  fi,  we  consider  the 
possible  location  of  fl  and  what  it  looks  like.  Note  that  for  every  k = 1, . . . , N 
the  following  representation  takes  place: 

J2n_  (j  - k)c  ■ 

x(ci,...,cN)  = k-  1+  3~^n = k-  1 + Coeffk  y(cu . . . ,cN),  (6) 

EJ=i  Ci 

where  CoeSk  = Y^Z  1 U ~ fy  Cj  + Yf=k+i  U ~ ty  ci  does  not  depend  on  ck. 
<0  >0 

If  only  Cfc  varies,  while  the  other  coefficients  are  fixed,  then  (x,  y)  lie  on  the 
straight  line  which  passes  through  the  point  ( k — 1, 0)  and  has  slope  1 /Coeffk- 
Here  CoeffMin  k < Coeffk  (ci, . . . , ck-\,ck+i,  ■■■,  cN)  < CoeffMaxk 
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Fig.  5.  Important  vertices  and  edges  in  the  space  of  parameters. 


for  every  choice  of  (ci, . . . ,Cfc_i,Cfc+i, . . . ,cw),  where 

CoeffMink  = Coeffk(cf , . . . , c£_vck+1, . . . ,c^),  . . 

CoefFMax k = Coefffc(q, . . . , ck_vc£+v  . . . , c£). 

The  domain  ft  lies  between  lines  with  maximal  and  minimal  slopes.  More 
precisely 

Lemma  4.  ft  lies  above  the  line  corresponding  to  CoefFMax  k iff 

CoefFMax k > 0 (otherwise  it  lies  below  the  line),  ft  lies  above  the  line  corre- 
sponding to  CoefFMink  iff  CoefFMin  k < 0 (see  Fig.  4). 

Let  us  now  introduce  the  following  notations  for  some  important  vertices 
and  edges  of  box  B.  There  are  2 N (among  2N)  important  vertices: 


M — (c+ , cj, . . . , c+_  j , ck  , ck+1, . . . , cN), 
P k — (CX  ! c2  I • • • I cfc-l>  cfc  > cit+l>  • • • 1 cyv)' 


(8) 


for  k = 1 ,N  and  2 N important  edges  connecting  these  vertices  (see  Fig.  5): 

ak{ck)  — (C1  i • ■ • ) ck-l’  Cfc>  cfc  + l ! . . . > civ  )'  ck  £ [cit^  > Ck  } ’ /g\ 

Pk{ck)  = (of , ■ ■ ■ , Ck_1,Ck,Ck+l,  . . . ,CN),  Ck  £ [cfc  , ck] . 

Edge  oik  connects  vertices  \k  and  A*+],  edge  f3k  - vertices  pk  and  Pk+i-  (Here 
An+i  = Pi  and  pn+i  — Ai). 

We  also  introduce  the  notation 


sum{a, . . . , cN)  — ci  + t-  cjv 


(10) 


which  will  be  useful  in  what  follows. 

As  explained  above,  it  is  clear  that  a*  is  mapped  onto  the  boundary  line 
of  ft  corresponding  to  CoefFMax  k,  and  (3k  is  mapped  onto  the  boundary  line 
corresponding  to  CoefFMin  k.  More  precisely,  we  have 

Lemma  5.  The  image  of  a*,  is  the  segment  between  fN(^k)  and  fiv(\k+i)  if 
y = 1 /sum(a(ck))  does  not  change  sign  while  Ck  € [ ck’ck  ]>  and  the  comple- 
ment of  the  straight  line  to  this  segment  otherwise.  In  other  words,  /,v  (ok ) 
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Fig.  6.  Image  of  a k depending  on  the  signs  of  y(Xk)  and  y(Xkyi). 


is  a segment  if  plane  c\  + C2  + • • • + cjv  = 0 does  not  intersect  edge  a*,  of  box 
B,  and  a complement  of  the  segment  otherwise  (see  Fig.  6). 

Of  course,  the  analogous  statement  holds  for  the  image  of  /?*..  In  what 
follows,  we  will  usually  formulate  only  statements  for  Afc  and  a*,,  and  omit  the 
analogous  ones  for  gk  and  /?*,. 

Conclusion  6.  /w(A *,)  and  /jv(^fc+i)  are  connected  by  segment  iff  they  lie 
in  the  same  (upper  or  lower)  half-plane,  i.e.  if  y(\k)  and  y(\k+i)  have  the 
same  signs. 

Note  further,  that  every  one  of  the  points  /yv(Afc)  (/yv(/ifc))  is  a point  of 
the  concatenation  of  two  boundary  segments  corresponding  to 
CoeffMaxk-i  and  CoeffMaxk  (CoeffMin k-i  and  CoeffMin k , respectively). 

In  order  to  make  this  precise  and  to  assure  that  all  boundary  can  be 
described  in  this  manner,  the  following  theorem  was  proved. 

Theorem  7.  /n(ci,  • ■ • , c/v)  belongs  to  the  boundary  of  domain  fl  iff  there 
exists  k = 1, . . . , N such  that  (ci, . . . , cN)  = ak{ck)  or  (ci, . . . , cN ) = /?fc(cfc) 
fOTCk  € [c)T,c+]. 

Again  the  proof  (which  is  relatively  long)  is  omitted  the  proof  uses  the 
“topological”  notion  of  point  neighborhood,  boundary  etc. 

Conclusion  8.  In  order  to  describe  the  boundary  of  the  domain  fi,  it  suffices 
to  move  along  the  following  contour  in  the  space  of  coefficients 


/ 

A2 

A3  — ► 

■ Ayv 

\ 

\ 

flN 

<—  PN-l  ■ ■ 

• ■*-  M2 

(11) 


Note  that  independently  of  the  specific  values  of  Cj  and  cj'  (j  = 1, ...  ,N), 
always  only  2N  (among  2N)  definite  vertices  and  edges  of  B and  in  definite 
order  participate  in  the  boundary  of  domain  ft. 


It  remains  to  “fill  in”  the  boundary  of  fl  with  51  itself.  Before  we  formulate 
the  main  result,  let  us  study  an  additional  property  of  domain  fl  and  its 
boundary. 
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Lemma  9.  The  domain  Q has  only  convex  boundary  vertices. 

Indeed,  it  was  shown  above  that  a boundary  vertex  of  fl  of  the  form 
/iv(^k)  is  the  intersection  of  the  lines  with  slopes  1 / CoeffMax  j.  __  j and 
1/ CoeffMax  k and  passing  through  points  ( k — 2, 0)  and  (k  — 1, 0)  respectively. 
Note  that  the  following  equation  holds: 

CoeffMaxk-i  — CoefFMaxk  = sum(Xk)  - 1/t/(A*)  (12) 

Assume  for  example  that  /^(Ak)  lies  in  upper  half-plane,  i.e.,  that  y(Xk)  > 0 
(the  analogous  consideration  can  be  done  for  the  lower  hyperplane).  Then 
CoeffMax k-i  > CoeffMax k , and  using  Lemma  4 we  get  that  in  any  one  of 
three  possible  cases  (see  Fig.  7)  this  vertex  is  convex. 

§5.  The  Main  Result  - Description  of  the  Domain  ft 

Theorem  10.  The  domain  ft  has  one  of  two  possible  forms  depending  on 
whether  sum( Ai)  = cf  + cj"  + • ■ • -fc)y  and  sum(pi)  = cj"  + + • • • + c^  has 

the  same  sign  or  not. 

Note  that  this  condition  is  equivalent  to  the  condition  that  the  plane 
Ci  + C2  + ■ • ■ + cjv  = 0 intersects  the  box  B in  the  space  of  coefficients. 

Case  l.  If  sum(Xi)  and  sum(pi)  have  the  same  sign,  then  fl  is  a convex 
bounded  polygon  inside  the  contour  (11),  (see  Fig.  8 and  9 - left  parts.  In  the 
figures  we  will  write  A*,  instead  of  /^(A*)  in  order  to  make  the  figures  clearer 
and  more  compact). 

Indeed,  let  us  assume  that  sum(pi)  > 0.  Then  y(Xk)  > 0 and  y(pk)  > 0 
for  every  k — 1, . . . , N (all  vertices  of  $1  lie  in  the  upper  half-plane).  According 
to  Conclusion  6,  the  boundary  of  fl  in  this  case  consists  of  segments  which 
form  a convex  (Lemma  9)  bounded  polygon. 

Case  2.  If  sum(pi)  < 0 and  sum{ X\)  > 0,  then  sum()  changes  its  sign  exactly 
once  when  upper  or  lower  chain  of  the  contour  (11)  is  traversed,  say  at  the 
segment  at  = [A* , A/t+i]  at  the  upper  chain  and  segment  fip  — [pp,  pp+i]  at 
the  lower  chain.  Then  the  domain  fi  is  a union  of  two  convex  unbounded 


Fig.  9.  An  example  of  the  domain  U in  the  3 D case. 

polygons.  The  first  is  inside  the  part  of  the  contour  that  belongs  to  the  upper 
half-plane,  and  the  second  is  inside  the  part  of  the  contour  that  belongs  to 
the  lower  half-plane. 

Upper  and  lower  parts  are  bounded  by  infinite  rays  with  the  same  slopes 
(corresponding  to  CoeffMaxk  and  CoeffMinp),  i.e.  it  is  a convex  bounded 
polygon  in  the  projective  plane  (see  Fig.  8 and  9 - right  parts). 


§6.  Example  of  an  Affine  Subspace  of  Lower  Dimension 

We  now  show  how  the  general  result  can  be  applied  to  the  construction  of 
p-flats  of  lower  dimensions,  for  the  “approximate”  line  in  3D  case.  For  ex- 
ample, if  we  would  like  to  describe  the  range  of  the  indexed  point  7fi3  = 

( 2c{13}  1 \ 

I -pffl-  jir,; , {i3}  {13)  I which  enters  in  the  representation  of  the  line,  then 

\C1  +C3  C1  +C3  ) 

we  can  reduce  it  to  the  consideration  of  fz  by  putting  C2  € [0, 0].  We  get  “void” 
connections  instead  of  edges  and  /?2,  and  finally  get  a convex  quadrilateral 
instead  of  convex  hexagon  (see  Fig.  10). 
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Fig.  10.  Range  of  W12  - the  first  indexed  point  corresponding  to  line  in  3D  case. 
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